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ABSTRACT 

Equations  for  the  determination  of  the  frequencies  of 
flexural  vibration  and  corresponding  modal  shapes  for  cylin- 
drically  aeolotropic  circular  plates  are  presented.  A  consid¬ 
erable  analysis  has  been  made  of  the  case  of  symmetric  vibra¬ 
tions  and  results  are  tabulated. 

The  theory  is  exemplified  by  an  analysis  of  the  vibration 
of  a  circular  plate  with  circular  stiffeners  for  which  the 
effective  elastic  compliances  had  been  determined  in  a  previous 
Investigation.  Nodal  patterns  and  corresponding  frequencies 
were  experimentally  determined  for  the  first  five  symmetric 
modes  and  the  first  four  asymmetric  modes  of  a  simply  supported 
stiffened  plate.  Comparison  between  the  theoretical  and  ex¬ 
perimental  results  for  the  first  three  symmetric  modes  is  pro¬ 


vided. 


Finally,  some  discussion  is  given  for  the  important 
problem  of  approximating  plates  with  attached  stiffeners  by 
equivalent  homogeneous  uniformly  thick  aeolotropic  plates. 

INTRODUCTION 

The  literature  contains  many  analyses  of  transverse  elas¬ 
tic  vibrations  of  plates  from  the  standpoint  of  small  deflec¬ 
tion,  thin-plate  theory1.  Past  research  deals  with  both  iso¬ 
tropic  and  orthotropic  plates  but  mainly  plates  which  are 

2 

rectangular  in  shape.  At  least  one  paper  in  the  recent  past 
deals  with  the  vibrations  of  circular  and  elliptical  ortho¬ 
tropic  plates  for  which  the  principal  axes  of  orthotropy  are 
parallel  or  perpendicular  to  a  given  diameter  in  the  circular 
plates  and  parallel  or  perpendicular  to  a  principal  diameter 
in  the  elliptical  plates. 

In  the  present  paper  both  a  theoretical  and  experimental 
study  is  made  of  the  flexural  vibration  of  circular  plates 
which  are  cylindrlcally  aeolotropic.  That  is  the  principal 
directions  of  aeolotropy  at  a  point  are  in  radial  and  trans¬ 
verse  directions.  Such  aeolotropy  may  occur  in  nature  as  in 

"^N.  J.  Huffington,  Jr.  and  W.  H.  Hoppmann  II,  "On  the  Trans¬ 
verse  Vibrations  of  Rectangular  Orthotropic  Plates",  J.  App. 
Mechanics,  25,  No.  3,  389-395  (1958). 

p 

W.  H.  Hoppmann  II,  "Flexural  Vibrations  of  Orthogonally 
Stiffened  Circular  and  Elliptical  Plates",  Proc.  of  3rd  U.S. 
Nat.  Cong,  of  App.  Mechanics,  181-187  (1958). 
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some  cuts  of  wood  or  may  be  manufactured  at  least  approxi¬ 
mately  in  the  case  of  plates  of  reinforced  materials  or  plates 
with  stiffeners  attached. 


DIFFERENTIAL  EQUATION  AND  ITS  SOLUTION 
FOR  THE  FREE  FLEXURAL  VIBRATION  OF 
CYLINDRI CALLY  AEOLOTROPIC  PLATES 

The  differential  equation  for  the  transverse  statical 

deflection  of  cylindrically  aeolotropic  plates  is  in  the 
'-i  4  c; 

literature  '  .  Its  extension  to  the  case  of  transverse 

flexural  vibration  is,  as  usual,  obtained  simply  by  adding  the 
inertia  term.  It  may  be  written  as  follows: 


D 


d  w 


dr 


De  d4w 

r^  d©'^ 


2D. 


d^w 


drd© 


D©  d2  w  +  2 
r^  dr2  r 


Q 


Dr©> 


d2w 


d© 


,  D©  dw  ,  7h  d2w  _  n 

+  — *■  H - yr  =  0 

r*3  dr  g  dt^ 


(1) 


^G.  F.  Carrier,  "The  Bending  of  the  Cylindrically  Aeolotropic 
Plate",  J.  App.  Mechanics,  vol.  11,  No.  3*  129  (19^*0  • 

^A.  M.  Sen  Gupta,  "Bending  of  a  Cylindrically  Aeolotropic  Cir¬ 
cular  Plate  with  Eccentric  Load",  J.  App.  Mechanics,  vol.  19* 
No.  1,  9  (1952). 

%.  G.  Lekhnitski,  "Anisotropic  Plates",  translated  by  E.  Z. 
Stowell,  American  Iron  and  Steel  Institute,  New  York,  June 
1956.  Contribution  to  the  Metallurgy  of  Steel,  No.  50,  p.22. 
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Equation  (l)  may  be  solved  by  the  method  of  separation 
of  space  and  time  variables.  Let  w  =  0(r,Q)elcot  and 
c-j2  ~  then  Eq.  (l)  becomes: 


2Dr  ^ 


-E  - 

r  Sr** 


S^0 

SrSO2 


r ^  Sr 


(2) 


where  a)  is  the  circular  frequency. 

Furthermore  we  may  likewise  separate  the  space  variables 
as  follows: 


0  = 


-  ^  Y(r), 


V=_«o 


Equation  (2)  then  becomes  an  ordinary  differential  equation 
on  Y(r)  as  follows: 

a4  2v\o  a2  .  v\  .  2Cr  d3  .  2v2Dre  a  Do  a2 

r  aT1  '  ~~7~  T”-  T  to3  ~7~  J7a?' 


ov2  DQ  j  2 

~T  <D9  +  Dr9>  +  4  i  -  <f  >  1  I(r) 
r  r  dr  Ct 


{d4  +  i.  d3  -  4  L?  +  Ij  D  ♦  %  -  a4}  Y(r) 

T»  T*  T*  Y*  * 
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where 


F  =  2v2o1  +  k2 


a4  =  -j£—  =  ,  and  G  =  v4k2  -  2V2(k2  +  o') 

°1  Dr  sDr 

The  latter  differential  equation  can  be  solved  in  series  form. 
Let 


Y(r) 


anr 


n+m 


(5) 


where  m  is  the  undetermined  index  of  the  series. 


Then  Eq.  (4)  gives  the  following  equation: 

[m(m-l)2(m-2)-m(m-2)F+G]a0rm+[  (l+m)  (m-l)m2-  (m-t-l)  (m-l)F+G]a1rm+1+ 

+  [(2+m)(l+m)2-m(2+m)F+G]a2rm+2 + 

+  [  (3+m)  (m+2  )2(m+l)- (m+1)  (S+nijF+Gja^r111^  + 

00 

+  ^  [  (n+m)(n+m-l)2(n+m-2)-F(n+m)(n+-m-2)+G]  anrn+m-a4an_^rn+m  =  0 

(6) 

Therefore,  if  Y(r)  is  to  be  an  integral  of  Eq.  (4),  the 
coefficient  of  each  power  of  r  in  Eq.  (6)  must  be  identically 
zero.  Setting  the  coefficient  of  rn+m  equal  to  zero,  the 
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recurrence  equation  for  the  coefficients  an  is; 

P 

[  (n+mXn+m-1)  (n+m-2  )-F(n+m)  (n+m-2  )+G]a  -a  a  k  =  0  n  >  4 

n  n-  i" 

(7) 

which  defines  all  the  coefficients  in  terms  of  four  of  them, 
say  ag  ,  ,  a2  ,  and  a^  as  soon  as  the  index  m  is 

known . 

If  aQ  0  ,  the  coefficient  of  rm  equated  to  zero 
gives  the  indicial  equation: 

m(m-l)2(m-2)  -  m(m-2)F  4-  G  =  0  (8) 


With  a  change  of  variable  m  =  x  +  1  ,  Eq.  (8)  can  be 
reduced  to  a  biquadratic  from  which  the  four  values  for  the 
indices  are: 


m  =  1 


(1+F)  +  J  (1+F)2  -  4(F+G)" 


(9) 


Once  the  elastic  compliances  of  a  certain  material  are 
known,  the  four  values  of  m  can  be  computed  from  Eq.  (9); 
the  recurrence  equation  for  every  root  can  then  be  obtained, 
and  the  series  solutions  for  that  particular  problem  obtained. 
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SOLUTIONS  FOR  THE  CASE  OF 
SYMMETRICAL  VIBRATIONS 

For  symmetrical  vibrations  the  differential  equation  is 

independent  of  0  and  all  terms  containing  v  vanish.  The 

o 

parameters  F  and  G  become  equal  to  k  and  0  respec¬ 
tively,  and  Eq.  (4)  reduces  to: 

[D4  +  JL  d3  _  ^  d2  +  £  D  _  a^]Y(r)  =  0  (10 ) 

r 

p 

Substituting  k  and  0  for  F  and  G  in  Eq.  (9), 
the  four  values  of  the  index  become: 


m.  =  0 

m2  =  2 

=  1  +  k 

=  1  -  k 

The  recurrence  equation  (7)  reduces  to: 


(ID 


(12) 


(13) 


The  coefficient  of  rm't‘1  ,  rm4’2  ,  require  that 

a-^  ,  a2  and  a^  be  equal  to  zero.  The  equation  correspond¬ 
ing  to  m  =  0  is: 


7 


a^a. 


a 


n 


n(n-2)[(n-l)2  -  k2] 


n  >  4 


and  gives; 


aYa, 


2*4(9-k  ) 


=  a*  =  a„  =  0 


6  “  d7 


a 


8 


a4a^ 


4  4 

a  *a  a. 


6-8(49-k  )  2-4.6-8(9-ir )(49-k2) 


and 


CiVr)  -  ,0[l  +  -JSSXL-  * - 


4  /  ,Q 

^  -  -V^*  L  4-  \CLT ) 

0L  2 • 4(9-k2 )  2-4-6-8(9-k2)(49-k2) 


+ 


& 
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2 . 4*6-  8*  10- 12  (9-k2 )  (49-k2 )  (121-k2 ) 


+ . + 


_  (ar)^ _  | 

22J'2jY(9-k^)(49-k^). .  •  [(4j-l)2-k2]  - 


] 


(14) 


Following  the  same  procedure  with  the  remaining  three 
Indices,  one  obtains; 
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02Y2(r)  -  a6p2[l  +  _ 


[ar)^  ,  (ar)® 


- i  - - -  + 

4*6(25-k  )  4 •  6 •  8 •  10 (2 5-k  ) (8l-k* ) 


+  - 5 + .  + 

4*  6  *8«  10*  12*  14(25-* )  (81-k^ )  (169-k  ) 


22^  (2j+l)!(25-k2)(8l-k2).  .  .  [(4j+l)2-k2] 


]  (15) 


CqYq(r)  =  a"r1+k[  1  +  - {SHE!! -  + 

33  0  L  8(2+k)(3+k)(5+k) 


+  _ j _  + 

8- l6(2+k) (3+k) (4+k) (5+k) (7+k) (9+k) 


+ - ^ -  +  . 

8-l6-24(2-t-k)(3+k)(4+k)(5+k)(6+k)(7+k)(9+k)(ll+k)(l3+k) 


+  ^i - - ]  (16) 

2JJ- y.  (2+k)(3+k). . .  (2j+k)(4j+l+k)(4j-l+k)  J 
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C4Y4(r)  =  a"'r1_k[  1  4-  - -  + 

44  0  L  8(2-k)(3-k)(5-k) 


l2£l 


8 


8.l6(2-k)(3-k)(4-k)(5-k)(7-k)(9-k) 


(ar) 


4J 


23J- j!(2-k)(3-k). . . (2j-k)(4j+l-k)(4j-l-k) 


1 


(17) 


By  using  the  ratio  test,  it  can  easily  be  shown  that  the 
four  solutions  above  converge  for  all  values  of  k  . 

The  functions  Y-^(r)  ,  Y2(r)  ,  Y^r)  and  Y4(r)  are 
solutions  corresponding  to  the  case  for  which  k  ,  the  square 
root  of  the  ratio  of  two  rigidities,  is  not  an  integer.  For 
the  case  k=l  ,  the  four  solutions  reduce  to  two. 

The  solution  is  then: 

w(r,  t)  =  [C1Y1(r>)  +  C2Y2(r)  +  C^r)  +  C4Y4(r )  ]ei£Ut  (l8) 


where  Y-^  ,  Yg  ,  Y^  ,  and  Y4  are  the  functions  obtained 
from  Eqs.  (l4)  through  (17)  and 


(19) 
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The  four  constants  0X  ,  C2  ,  ,  and  can  be 

determined  from  the  boundary  conditions  for  an  annulus.  For 
a  solid  plate  Y^  is  inadmissible  because  of  its  singularity 
at  r  =  0.  For  this  case  we  write  =  0.  Also  for  sym¬ 
metrical  vibrations  it  can  be  shown  that  C2  =  0.  The  shear 
in  the  plate  at  radius  r  must  balance  the  inertia  forces  as 
follows: 


2ttt  Qr 


a2w 

S7 


dS 


(20) 


where 

are: 

w(r,t) 

d2w 


Qr  is  the  shear  and  the  displacement  and  acceleration 

=  (G1Y1  4-  C2Y2  +  C3Y3)ela>t 

-  a)2(C1Y1  +  C2Y2  +  C3Y3)ela)t 


The  shear  is  given  by: 


%  -  -  U>r( 


d^w  +  _1_  d2w  \ 


r'  a? 


r  dr' 


©  dw  i  io>t 
—  Je 

r  dr 


and  so  we  have  as  a  condition  to  be  satisfied 


d^w 


d2w 


2 

uw  ,  uw  k  dw  i 

-  L  r  — -  —  —  J 

drJ  dr  r  dr 


(21) 
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Substituting  the  values  of  the  derivatives  of  the  first  four 
terms  of  Y^  ,  Y2  and  Y^  into  Eq.  (21),  we  obtain  finally: 

-  C2[2(l^k2)]  =  0 

and  C2  =  0  ,  since  k  1  .  Hence,  the  solution  for  a  solid 
plate  becomes: 

w(r,t)  =  [C1Y1  +  C3Y3]ela5t  (22) 


where  the  two  constants,  and  C3  ,  will  have  to  be  deter¬ 

mined  from  the  remaining  boundary  conditions. 

1.  Clamped  Plate 

The  boundary  conditions  for  a  plate  clamped  at  r  =  a 

are: 

w  =  0  at  r  =  a 

—  =  0  at  r  =  a 

dr 


Inserting  these  conditions  iryto  Eq.  (22),  we  obtain: 


C1Y1(aa)  +'C3Y3(cxa)  =  0 

C^^aa)  +  C3Y3(aa)  =  0 


(23) 


where  the  prime  indicates  the  first  derivative  with  respect 
to  r  .  Eliminating  C1  and  C3  from  Eq.  (23)  above  leads 
to  the  frequency  equation: 


12 


0 


(24) 


Y1(aa)Y^(oa)  -  Y3(aa)Y^(cxa)  = 

The  roots  of  Eq.  (24)  give  the  values  of  a  ,  which  when 
inserted  in  Eq.  (19)  yield  the  values  of  the  frequencies  of 
the  clamped  plate. 

2 .  Simply  Supported  Plate 

For  the  simply  supported  plate  the  boundary  conditions 


become : 

w  =  0 

at  r  =  a 

,2 

M  - 

+  _£*'=  0  at  r  =  a 

r  dr2 

r  dr 

Using 

the  boundary  conditions  the  frequency  equation 

becomes: 


v  V 

Y-j^aa)  [Y3(aa)  +  —  Y3(aa)]  -  Y^(aa )  [Y^(aa)  -f~Y|(aa)]  =  0 

3.  at 

(25) 

Again,  the  roots  of  the  frequency  equation  (25)  are  to 
be  inserted  into  Eq.  (19)  to  determine  the  frequencies  of  the 
simply  supported  plate. 


The  first  six  terms  of  the  infinite  series  were  used  in 
the  expansion  of  Eq.  (25) >  and  the  first  three  frequencies 
were  calculated  for  several  sets  of  assumed  values  of  the 
elastic  constants.  The  results  are  shown  in  Table  I. 
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EXPERIMENTS  ON  CIRCULARLY  STIFFENED 
CIRCULAR  PLATES 


The  apparatus  for  determining  the  frequencies  of  vibration 
and  the  nodal  patterns  is  shown  in  Figure  1.  The  plate  under 
study  was  driven  at  resonance  by  an  electromagnet  mounted  be¬ 
neath  it  on  a  solid  test  stand.  A  rectangular  piece  of  lami¬ 
nated  iron  was  used  in  order  to  develop  a  driving  force.  It 
was  cemented  to  the  plate  through  a  flat  circular  washer  in 
order  to  create  a  symmetrical  point  of  dynamic  loading. 

The  frequencies  and  nodal  patterns  were  determined  by 
means  of  a  capacity  type  of  pickup  which  can  be  moved  across 
the  surface  of  the  plate.  A  nodal  line  is  detected  by  the 
change  in  shape  of  a  Lissajous  figure  on  a  cathode-ray  oscillo¬ 
scope  to  which  the  electrical  output  from  the  pickup  is  fed 
through  a  preamplifier.  On  account  of  the  pull-pull  nature 
of  the  magnet,  the  frequency  from  the  oscillator  is  one  half 
the  frequency  of  the  driving  force.  The  output  from  the  os¬ 
cillator  Is  connected  to  one  set  of  plates  of  the  cathode-ray 
oscilloscope,  and  the  magnetic  drive  is  connected  across  the 
other  set  of  plates. 

Since  the  frequency  ratio  is  2  to  1,  a  figure  eight 
Lissajous  is  shown  on  the  screen  of  the  scope.  Crossing  the 
nodal  line  with  the  probe  or  pickup  turns  the  figure  eight 
Lissajous  into  a  horseshoe-shaped  figure  on  one  side  of  the 
line,  and  inverts  the  figure  on  the  opposite  side.  This  fact 
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provides  a  precise  means  of  constructing  the  nodal  pattern. 

The  experimental  frequencies  and  nodal  patterns  for  the 
plate  under  study  are  shown  in  Figure  2.  The  theoretical  and 
experimental  frequencies  for  the  first  three  symmetrical  modes 
of  vibration  are  shown  in  Table  II.  The  theoretical  frequen¬ 
cies  were  obtained  from  Table  I  for  the  elastic  constants 
determined  from  previous  static  tests. 

Table  II 

Comparison  of  Experimentally  Determined 
and  Calculated  Frequencies  for 
Symmetrical  Vibration 


h 

c .p. s. 

f2 

c.p. s. 

f3 

c .p. s. 

Experimental 

460 

1634 

3830 

Theoretical 

392 

1990 

4708 

While  the  discrepancy  between  the  theoretical  and  experi¬ 
mental  results  is  fairly  large,  it  may  be  noted  that  the 
theoretical  frequencies  were  based  on  experimentally  deter¬ 
mined  elastic  compliances  for  the  stiffened  plate  material 
and  therefore  subject  to  some  error.  Also,  the  rotatory 
moment  of  inertia  is  not  included  in  the  theoretical  solution. 
Its  effect  is  to  reduce  the  higher. frequencies  as  pointed  out 
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in  a  paper  by  Thorkildsen  and  Hoppmann  . 

A  much  better  agreement  with  experimental  results  is  to 
be  expected  if  the  stiffeners  are  more  closely  spaced  and 
symmetrically  disposed  with  respect  to  the  center  plane  of  the 
plate,  thus  providing  a  better  approximation  to  the  condition 
of  homogeneity  assumed  in  the  theory. 


DISCUSSION  AND  CONCLUSIONS 

The  reason  for  the  discrepancy  between  theory  and  experi¬ 
ment  in  the  present  study  can  reasonably  be  attributed  to  the 
inappropriateness  of  the  design  of  the  experimental  plate 
which  does  not  sufficiently  well  meet  the  requirement  of  homo¬ 
geneity  on  which  the  theory  is  based.  It  is  clear  that  homo¬ 
geneity  will  be  better  approximated  if  a  larger  number  of 
appropriately  proportioned  stiffeners  are  used.  Since  the 
radius  of  the  stiffener  decreases  for  stiffeners  closer  to  the 
center  of  the  plate,  the  cross-sections  must  be  likewise  re¬ 
duced.  The  cross-section  of  the  stiffeners  used  in  the  pres¬ 
ent  investigation  were  reduced  in  accord  with  a  rational  plan 
but  it  is  considered  that  future  research  is  necessary  to 
answer  this  question  in  a  satisfactory  manner. 

R.  L.  Thorkildsen  and  W.  H.  Hoppmann  II,  "Effect  of  Rotatory 
Inertia  on  the  Frequencies  of  Vibration  of  Stiffened  Plates", 
J.  Appl.  Mechanics,  298  ( 1959 ) • 


17 


It  is  also  clear  that  if  the  stiffeners  were  formed 
symmetrically  with  respect  to  the  middle  surface  and  rotatory 
Inertia  taken  into  account  the  agreement  between  the  theo¬ 
retical  and  experimental  results  would  definitely  be  improved. 
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